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Abstract. We prove that many aspects of the differential geometry of embed- 
ded Riemannian manifolds can be formulated in terms of multi linear algebraic 
structures on the space of smooth functions. In particular, we find algebraic 
expressions for Weingarten's formula, the Ricci curvature and the Codazzi- 
Mainardi equations. 

For matrix analogues of embedded surfaces we define discrete curvatures 
and Euler characteristics, and a non-commutative Gauss-Bonnet theorem is 
shown to follow. We derive simple expressions for the discrete Gauss curva- 
ture in terms of matrices representing the embedding coordinates, and a large 
class of explicit examples is provided. Furthermore, we illustrate the fact that 
techniques from differential geometry can carry over to matrix analogues by 
proving that a bound on the discrete Gauss curvature implies a bound on the 
eigenvalues of the discrete Laplace operator. 
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1. Introduction 

It is generally interesting to study in what ways information about the geometry of 
a differentiable manifold £ can be extracted as algebraic properties of the algebra 
of smooth functions C°°(E). In case £ is a Poisson manifold, this algebra has 
a second (apart from the commutative multiplication of functions) bilinear (non- 
associative) algebra structure, the Poisson bracket. The bracket is compatible with 

l 
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the commutative multiplication via Leibniz rule, thus carrying the basic properties 
of a derivation. 

On a surface E, with local coordinates u 1 and u 2 , one can define 

' \f9\du 1 du 2 du 1 du 2 

where g is the determinant of the induced metric tensor, and one readily checks that 
(C X (E), {•,•}) is a Poisson algebra. Having only this very particular combination 
of derivatives at hand, it seems at first unlikely that one can encode geometric 
information of E in Poisson algebraic expressions. Surprisingly, it turns out that 
many differential geometric quantities can be computed in a completely algebraic 
way, cp. Theorem 13 . 71 and Theorem 13.171 For instance, the Gaussian curvature of 
a surface embedded in R m can be written as 



(1.1) K= ^ (^{{x j ,x k },x k }{{xi,x l },x l }--{{x>,x«},x l }{{x>,x%x 1 } 
j,M=i V 

where X % {u , u ) are the embedding coordinates of the surface. 

For a general n-dimensional manifold E, we are led to consider Nambu brackets 
|Nam73] . i.e. multi-linear alternating n-ary maps from C' X (E) x •■• x C X (E) to 
C X (E), defined by 

{/l, ■•-,/»} = ^£ ai - Q " (dajl)---(dajn). 

\J9 

In the case of surfaces, our initial motivation for studying the problem came from 
matrix regularizations of Membrane Theory. Classical solutions in Membrane The- 
ory are 3-manifolds with vanishing mean curvature in R 1 '^. Considering one of the 
coordinates to be time, the problem can also be formulated in a dynamical way 
as surfaces sweeping out volumes of vanishing mean curvature. In this context, a 
regularization was introduced replacing the infinite dimensional function algebra 
on the surface by an algebra of N x N matrices [GH82 . If we let T a be a linear 
map from smooth functions to hermitian N a x N a matrices, the main properties of 
the regularization are 



lim \\T a (f)T a (g)-T a (fg)\\=0, 



lim 



rj—[T a (f), T a (h)] — T a ({f, h}) 



0, 



where H a is a real valued function tending to zero as N a — > 00 (see Section @] for 
details), and therefore it is natural to regularize the system by replacing (commu- 
tative) multiplication of functions by (non-commutative) multiplication of matrices 
and Poisson brackets of functions by commutators of matrices. 

Although we may very well consider T a (j£;), its relation to T a (f) is in general 
not simple. However, the particular combination of derivatives in T a ({f, h}) is ex- 
pressed in terms of a commutator of T a (f) and T a (h). In the context of Membrane 
Theory, it is desirable to have geometrical quantities in a form that can easily be 
regularized, which is the case for any expression constructed out of multiplications 
and Poisson brackets. For instance, solving the equations of motion for the reg- 
ularized membrane gives sequences of matrices that correspond to the embedding 
coordinates of the surface. Since the set of solutions contains regularizations of 
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surfaces of arbitrary topology, one would like to be able to compute the genus cor- 
responding to particular solutions. The regularized form of provides a way of 
resolving this problem. 

The paper is organized as follows: In Section[2]we introduce the relevant notation 
by recalling some basic facts about submanifolds. In Section [3] we formulate several 
basic differential geometric objects in terms of Nambu brackets, and in Section [3. II 
we provide a construction of a set of orthonormal basis vectors of the normal space. 
Section 13.21 is devoted to the study of the Codazzi-Mainardi equations and how 
one can rewrite them in terms of Nambu brackets. In Section 13.41 we study the 
particular case of surfaces, for which many of the introduced formulas and concepts 
are particularly nice and in which case one can construct the complex structure in 
terms of Poisson brackets. 

In the second part of the paper, starting with Section[4j we study the implications 
of our results for matrix regularizations of compact surfaces. In particular, a discrete 
version of the Gauss-Bonnet theorem is derived in Section 14.11 and a proof that the 
discrete Gauss curvature bounds the eigenvalues of the discrete Laplacian is found 
in Section 14.41 

2. Preliminaries 

To introduce the relevant notations, we shall recall some basic facts about submani- 
folds, in particular Gauss' and Weingarten's equations (see e.g. [KN96a. KN96b for 
details). For n ^ 2, let E be a n-dimensional manifold embedded in a Riemannian 
manifold M with dim M = n + p = m. Local coordinates on M will be denoted by 
x 1 , . . . ,x m , local coordinates on S by u 1 , . . . , and we regard x 1 , . . . , x m as being 
functions of tt 1 , . . . , u n providing the embedding of E in M. The metric tensor on 
M is denoted by g^ and the induced metric on E by g a b', indices i,j,k,l,n run from 
1 to m, indices a,b,c,d,p,q run from 1 to n and indices A,B,C,D run from 1 to 
p. Furthermore, the covariant derivative and the Christoffel symbols in M will be 
denoted by V and Tj k respectively. 

The tangent space TE is regarded as a subspace of the tangent space TM and 
at each point of E one can choose e a = (d a x l )di as basis vectors in TE, and in this 
basis we define g a b = <?(e a ,e&)- Moreover, we choose a set of normal vectors Na, 
for A = l,...,p, such that §(Na, Nb) = Sab and §(Na, e a ) = 0. 

The formulas of Gauss and Weingarten split the covariant derivative in M into 
tangential and normal components as 



where X, Y £ TE and V X Y, W A (X) 6 TE and a(X,Y), D X N A £ TE- 1 . By 
expanding a(X, Y) in the basis {N±, . . . , N p } one can write (|2.ip as 



and we set /ia,o& = ^A(e a , e&)- From the above equations one derives the relation 



(2.1) 
(2.2) 



S7 X Y = V X Y + a(X,Y) 
V X N A = -W A (X)+D X N A 



v 



(2.3) 




(2.4) 



hA,ab = -g(e a ,V b N A ), 
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as well as Weingarten's equation 

(2.5) h A {X,Y)=g(W A (X),Y), 

which implies that (W,i)£ = g ac h AiC b, where g ab denotes the inverse of g ab - 

From formulas (|2.1[) and (|2.2j) one obtains Gauss' equation, i.e. an expression 
for the curvature R of E in terms of the curvature R of M, as 

g(R(X,Y)Z,V) = g(R(X,Y)Z,V) - g(a(X,Z),a(Y,V)) 
+ g(a(Y,Z),a(X,V)), 

where X,Y,Z,V e TE. As we shall later on consider the Ricci curvature, let us 
note that (|2.6j) implies 

(2.7) Ul = gP d g ac g(R(e C: e d )e b , e a ) + £ [(W A ) a a (W A ) p b - (W 2 A f b ] 

A=l 

where 1Z is the Ricci curvature of E considered as a map TE — » TE. We also recall 
the mean curvature vector, defined as 

1 p 

(2.8) H=-V (tr W A )N A . 



3. Nambu bracket formulation 

In this section we will prove that one can express many aspects of the differential 
geometry of an embedded manifold E in terms of a Nambu bracket introduced on 
C°°(E). Let p : E — » M. be an arbitrary non- vanishing density and define 

(3-1) {/!,...,/»} = - p S ai - a "(daJl)---(daJn) 

for all fx,.-. , f n 6 C X (E), where e ai "' a ' n is the totally antisymmetric Levi-Civita 
symbol with e 12 "' n = 1. Together with this multi-linear map, E is a Nambu-Poisson 
manifold. 

The above Nambu bracket arises from the choice of a volume form on E. Namely, 
let ui be a volume form and define {/i, . . . , /„} via the formula 

(3.2) {/i,...,/„}cj = dfx a ••• a df n . 

Writing u> = pdu 1 a ■ • ■ a du n in local coordinates, and evaluating both sides of 
(|3.2p on the tangent vectors d u i , . . . , d u ™ gives 

■ ■ ■ ■ « - \ *' (fc^) - V"" ^ ■ ■ ■ 

To define the objects which we will consider, it is convenient to introduce some 
notation. Let a; 1 (u 1 , . . . , u n ), . . . , x m (u 1 , . . . , u n ) be the embedding coordinates of 
E into M, and let n^(w, 1 , . . . ,u n ) denote the components of the orthonormal vectors 
N Al normal to TE. Using multi-indices I = i\ ■ ■ ■ i„_i and a = a± ■ ■ ■ a„_i we define 

{/,*'} sttzS^...,**- 1 } 
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together with 

dsx 1 = (d ai x^) (d a2 x 12 ) ■ ■ • (5 a „_ 1 ^- 1 ) 

iysfiA) 1 = (V ai N A ) 11 (V a2 N A y 2 ■ ■ ■ (^ ttn _ 1 N A ) in " 1 

9l.J = 9i 1 j 1 9i 2 j 2 ' ' ' 9i n _ 1 j n _ 1 
9ac = 9a 1 c 1 9a 2 c 2 ' ' ' 3a,,-ic„_i • 

We now introduce the main objects of our study 

1 1 F aS 

^ s * = J {n l ly — {daxl) 

(3.5) r i J = ^=TyiV ( ^ /)( ^ )J 

from which we construct 

(3.6) (P 2 Y k = V u V kJ g u 
(3-7) (B A ) lk = V u {TA) Jk - 9lJ 
(3-8) (S A T A y k = (S A f(T A ) Jk g IJ . 

By lowering the second index with the metric g, we will also consider V 2 , B A and 
T A S A as maps TM — » TM. Note that both S A and T A can be written in terms of 
Nambu brackets, e.g. 



T Ti_ ("!)' 

'A 



Let us now investigate some properties of the maps defined above. As it will appear 
frequently, we define 

(3.9) 7 =— ■ 

P 

ft is useful to note that (cp. Proposition 13. 3p 

7 2 = Yi -9 l3 { x \x I }g I j{x\x J }, 
■■j.i.i i "• 

and to recall the cofactor expansion of the inverse of a matrix: 
Lemma 3.1. Let g ab denote the inverse of g ab and g = det(g a b). Then 

(3.10) gg ba = —^^■■■^- 1 ^---W-x g g . . . 9an _ ihn _ x . 

(n — ly. 

Proposition 3.2. For X e TM it holds that 
(3.H) V 2 (X)= 1 2 g(X,e a )g ab e b 

(3.12) B A {X) = -l 2 g(X, V a N A )g ah e b 

(3.13) S A T A (X) = 7 2 (det W A )g(X, V a N A )hfe b , 
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and for Y e TE one obtains 

(3.14) V 2 (Y)=j 2 Y 

(3.15) B A (Y) = l 2 W A {Y) 

(3.16) S a Ta{Y) = - 7 2 (det W A )Y. 

Proof. Let us provide a proof for equations (|3.1ip and f|3 . 14|) : the other formulas 
can be proved analogously. 

V\X) = ■p u V iJ g IJ g jk X k d i = _ fax 1 ) fax 1 ) {d c x>) (d s x J )g IJ g jk X k d i 

9 \ n *-)• 



ySaiC! ■ ■ ■ 5a„_ lC „_! (<9 a x l ) (d c x 3 )g jk X d. 



p^(n — 1)! v 7 v y J 

= 7 V c (^ i ) {dc^)g jk X k di = l 2 g(X, e c )g ca e a . 

Choosing a tangent vector Y = Y c e c gives immediately that T 2 (Y) = j 2 Y. □ 

For a map B : TM — > TM we denote the trace by Tr£> = B\ and for a map 
W : TE — > TE we denote the trace by tr W = W£. 

Proposition 3.3. It holds that 

(3.17) lTrP 2 = 7 2 

n 

(3.18) Tr B A = I 2 trW A 

(3.19) - TtSaTa = -7 2 (detT^A). 

n 

Remark 3.4. For a hypersurface (with normal N = n l di) in R™ +1 , 

n 1*^ 1 ■, - • - ■, X n }{rii 1 , . . . , ?7^ n } 



(3.20) detW = (-l) 



>,:>■" 



}{x kl , ■ ■ -,x kn 



7«! 

the signed ratio of infinitesimal volumes swept out on S n (by N), resp E (which 
can easily be obtained directly by simply writing out the determinant of the second 
fundamental form, h = det(— d a x l d^rii)); in fact, all the symmetric functions of the 
principal curvatures are related to ratios of products of two Nambu brackets (cp. 
the paragraph after Proposition I3.1lj) . Namely, the /c'th symmetric curvature is 
given by 



(3.21) (-1) 



x 



, x n , • ■ • , rii k , Xi k+1 , • ■ • , Xi n } 



{x kl , . . . ,x k "}{x kl , ■ ■ ■ ,x kn } 



A direct consequence of Propositions l3T21 and l373l is that one can write the projection 
onto TE, as well as the mean curvature vector, in terms of Nambu brackets. 

Proposition 3.5. The map 

(3.22) 7- 2 ^ 2 = tt^V 2 : TM -> TE 

TrV 
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is the orthogonal projection of TM onto TE. Furthermore, the mean curvature 
vector can be written as 

1 A = l 

Proposition l3.2l tells us that j^ 2 B A equals the Weingarten map Wa, when restricted 
to TE. What is the geometrical meaning of B A acting on a normal vector? It turns 
out that the maps Ba also provide information about the covariant derivative in 
the normal space. If one defines (Dx)ab through 

v 

D x N A =Y i (Dx)abN b 

B = l 

for X e TE, then one can prove the following relation to the maps Ba- 
Proposition 3.6. For X e TE it holds that 

(3.23) g(B B {N A ),X)= 1 2 (D x ) AB . 

Proof. For a vector X = X a e a , it follows from Weingarten's formula (12.2[) that 

(Dx)ab = g(V x N A ,N B ). 

On the other hand, with the formula from Proposition 13.21 one computes 

g(B B (N A ).X) = - 1 2 g(N A ,V a N B )g ab g bc X c = -j 2 g(N A ,y x N B ) 

= - 1 2 {D X ) BA =1 2 {D x ) ab . 

The last equality is due to the fact that D is a covariant derivative, which implies 
that = D x g(N A ,N B ) = g(D x N A , N B ) + g(N A , D X N B ). □ 

Thus, one can write Weingarten's formula as 

p 

(3.24) l 2 V x N A = -Ba(X) + J] g(B B (N A ),X)N B , 

B = l 

and since Ha(X, Y) = 7 -2 <7(£>a(-'0, ^0 Gauss' formula becomes 

1 p 

(3.25) V X Y = V X Y + — g(B A (X),Y)N A . 

1 A=l 

Let us now turn our attention to the curvature of E. Since Nambu brackets involve 
sums over all vectors in the basis of TE, one can not expect to find expressions 
for quantities that involve a choice of tangent plane, e.g. the sectional curvature 
(unless E is a surface). However, it turns out that one can write the Ricci curvature 
as an expression involving Nambu brackets. 

Theorem 3.7. Let 1Z be the Ricci curvature o/S, considered as a map TE — > TE ; 
and let R denote the scalar curvature. For any X e TE it holds that 

(3.26) R{X) = ^(r 2 y k (p 2 ) lm R l3kl X^d m + UTrB A )B A (X) - B\{X)] 

' > A=l 

(3.27) R = \(V 2 f(V 2 fR l]kl + [(TtBa) 2 - TrB 2 A (X) , 

' ' A=l 

where R is the curvature tensor of M. 
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Proof. The Ricci curvature of E is defined as 

K=9 aC 9 Pd 9(R(ec,e d )e bl e a ) 
and from Gauss' equation (|2.6j) it follows that 



K = g pd g ac g(R(e c ,e d )e b ,e a )+g ac gP d j] (h AM h A<ac - h A>bc h AM \ 

A = l 

Since (T^a)^ = g ac h A ,cb one obtains 

K = g ac g pd g(R(e c ,e d )e b ,e a ) + j] \(tvW A )(W A f b - {W 2 A fX 

A = l 

and as B A {X) = -/ 2 W A (X) for any X e TT,, and Tr Ba = 7 2 trW A , one has 
K(X) = g ac gP d g(R(e c ,e d )e b ,e a )X»e p + \ j] [( Tr B A )B A (X) - B\{X) 

^ A = l 

By expanding the first term as 

g ac gP d X b R ijk l (dax 1 ) {d b x j ) (8 c x k ) (d d x l ) (d p x m )d m 

= g2{n \ 1)}2 s^e dd g 6 aS e^g ss X b R ijkl (d a x l ) (d b x>) (d c x k ) (d d x l ) (d p x m )d m 

= ... = ^{V 2 ) lk {V 2 ) lm R ljkl X^ m 
one obtains the desired result. □ 

3.1. Construction of normal vectors. The results in Section [3] involve Nambu 
brackets of the embedding coordinates and the components of the normal vectors. 
In this section we will prove that one can replace sums over normal vectors by sums 
of Nambu brackets of the embedding coordinates, thus providing expressions that 
do not involve normal vectors. 

It will be convenient to introduce yet another multi-index; namely, we let a = 
i\ . . . i v -\ consist of p — 1 indices all taking values between 1 and m. 

Proposition 3.8. For any value of the multi-index a, the vector 

(3.28) Z a = - 1 =- g i ^ jkl ... kna {x k \...,x k -}8 i , 

7(nY(p-l)!J 

where £i x —i m is the Levi-Civita tensor of M, is normal to TT,, i.e. g(Z a ,e a ) = 
for a = 1,2, ... ,n. For hyper surf aces (p = 1), equation iS. 28\) defines a unique 
normal vector of unit length. 

Proof. To prove that Z a are normal vectors, one simply notes that 

1 (n\^(p-\y)g{Z a ,e a ) = jU" 1 " "^...^ (d a x?) (d ai x kl ) ■ ■ ■ (d an x k ") = 0, 
since the n + 1 indices a, a\, . . . , a n can only take on n different values and since 



{d a x^){d ai x kl ) ■ ■ ■ (da n x kn ) is contracted with ejk 1 .--k n a which is completely anti- 



symmetric in j, ki,..., k n . Let us now calculate \Z\ = g(Z, Z) when p = 1. Using 



thatQ 



one obtains 



\ Z \ 2 = -Tvi5hi> ■ ■ ■ 9i n v e ikl ... k J h - 1 - {x k - , . . . , x fe " Kz'i , . . . , a: 1 - 



= 7^{ xll '---' xln }9hli-'-9u'Jx l \...,x l 'n} 

= -^(n-iy.TrT 2 = -^-(n-l)\n 7 2 = 1, 
which proves that Z has unit length. □ 

If the codimension is greater than one, Z a defines more than p non-zero normal 
vectors that do not in general fulfill any orthonormality conditions. In principle, 
one can now apply the Gram-Schmidt orthonormalization procedure to obtain a set 
of p orthonormal vectors. However, it turns out that one can use Z a to construct 
another set of normal vectors, avoiding explicit use of the Gram-Schmidt procedure; 
namely, introduce 

Zi=g{Z a ,Z% 

and consider it as a matrix over multi-indices a and /3. As such, the matrix is 
symmetric (with respect to g a p = §i 1 j 1 • • ■ g~i _ t j _J and we let E a ^ , /i Q denote or- 

J P 

Using these eigenvectors to define 

N a = EP a Zp 



thonormal eigenvectors (i.e. gg a E^E^ = 8 a p) and their corresponding eigenvalues. 



one finds that g{N a ,Nf}) = /J, a S a p, i.e. the vectors are orthogonal. 
Proposition 3.9. For Z? = g ZJ Z l a Z jl3 it holds that 

(3.29) Z s a Zf = Z£ 

(3.30) Z a a = p. 

Proof. Both statements can be easily proved once one has the following result 

(3.31) Z l a Z> a =g lJ - ^(V 2 Y\ 
which is obtained by using that 

E*fc 1 ...*„ae ,,1 " J » a = (p-l)l(4i^---^). 

Formula p.30[) is now immediate, and to obtain Q3.29P one notes that since Z a e 
TE X it holds that V 2 (Z a ) = 0, due to the fact that V 2 is proportional to the 
projection onto TS. □ 



4n our convention, no combinatorial factor is included in the anti-symmetrization; for instance, 
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From Proposition 13.91 it follows that an eigenvalue of Z is either or 1, which 
implies that N a = or g(N a , N a ) = 1: and that the number of non-zero vectors is 
Tr Z = Z" = p. Hence, the p non-zero vectors among N a constitute an orthonormal 
basis of TT,- 1 , and it follows that one can replace any sum over normal vectors Na 
by a sum over the multi-index of N a . As an example, let us work out some explicit 
expressions in the case when M = M. m . 

Proposition 3.10. Assume that M = W n and that all repeated indices are summed 
over. For any X e TE one has 



(3.32) J (TrB A )B A (Xy = ——-W k {{^, x J }, x J }{x\ x ! }{X\ x 1 } 

(3.33) j] B\(Xy = 1 IP k {x\ x x J }{X k , S J ), x 1 } 

'(»-')! 



(—11™ 

(3.34) ^ ( Tr Ba) N\ = A—Uw k {{x\ x 1 }, x 1 



wh 



ere 



1 '— 2\V 



(3.35) IF = S* - -j (P 2 ) 

is the projection onto the normal space. 

Proof. Let us prove formula (|3.32[) ; the other formulas can be proven analogously. 
One rewrites 

(TrB A )B A (Xy = ^ {x^,x J }{x J ,n j A }{x i ,x I }{x I ,n k A }X k 

= {n ^ l)[2 n\n k A {x J 1 {x^x J }}{x\x I }{x I 1 X k } 

since n J A {x : ',x J } = n A X k = 0, due to the fact that Na is a normal vector. By 
replacing n A n k A with N^N k and using the fact that 



A A 

1 



2\U 



KK = s 13 - -j (p 2 ) 

one obtains 

(TrB A )B A (Xy = * n> k {{xt, x J }, x J }{x\ z J }pf fe , x 1 }. □ 

For hypersurfaces in K™ +1 , the "Theorema Egrcgium" states that the determinant 
of the Weingarten map, i.e the "Gaussian curvature" , is an invariant (up to a 
sign when E is odd-dimensional) under isometries (this is in fact also true for 
hypersurfaces in a manifold of constant sectional curvature). From Proposition 13.31 
we know that one can express det W a in terms of Tr SaTa- 

Proposition 3.11. Let E be a hypersurface in M. n+1 and let W denote the Wein- 
garten map with respect to the unit normal 

Z = —g l 3 £jkK {x k 1 X K }. 

7n! 
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Then one can write det W as 
1 

7(771!) 



x {z\ , x Kl }, . . . , {x*- 1 , x^}}{x J ,{x l ,x^}}. 

In fact, one can express all the elementary symmetric functions of the principle cur- 
vatures in terms of Nambu brackets as follows: The elementary symmetric functions 
of the eigenvalues of W is given (up to a sign) as the coefficients of the polynomial 
det(W-tl). Since B(X) = for all X e TE- 1 and B(X) = j 2 W(X) for all X e TE, 
it holds that 

-tdet(W-tt n ) =det( 7 - 2 e-rt„ +1 ) = det(g-^ 7 2 l» +1 ) 

which implies that the coefficient of t k in det(IT — £1) is given by the coefficient of 
t k+i in -det(B-i7 2 l) 7 2 (™-' £ ). 

3.2. The Codazzi-Mainardi equations. When studying the geometry of em- 
bedded manifolds, the Codazzi-Mainardi equations are very useful. In this section 
we reformulate these equations in terms of Nambu brackets. 

The Codazzi-Mainardi equations express the normal component of R(X, Y)Z in 
terms of the second fundamental forms; namely 

g (R{X, Y)Z, N A ) = {V x h A ){Y,Z) - (S7 Y h A )(X,Z) 



(3.36) 



p 



J \g{D x N B ,N A )h B {Y,Z) - g{D Y N B ,N A )h B {X,Z)\ 

A = l 

for X,Y,Ze TE and A = I,..., p. Defining 

W A {X,Y) = (VxW A )(Y) - (V y Wa)(X) 



(3.37) 



p 



J \g(D x N B ,N A )W B (Y) - g(D Y N B ,N A )W B (X)] 

B = l 

one can rewrite the Codazzi-Mainardi equations as follows. 

Proposition 3.12. Let II denote the projection onto TE X . Then the Codazzi- 
Mainardi equations are equivalent to 

(3.38) Wa(X, Y) = -(1 - IT) (R(X, Y)N A ) 
forX,YeTZ and A = l,...,p. 

Proof. Since h A (X, Y) = g(W A (X), Y) (by Weingarten's equation) one can rewrite 
(j3~36l) as 

(3.39) g(W A (X,Y),Z) = g(R(X,Y)Z, N A ) , 
and since g(R(X,Y)Z,N A ) = -g(R(X,Y)N A , Z) this becomes 

(3.40) g(W A (X, Y) + R(X, Y)N A , Z) = 0. 
That this holds for all Z e TE is equivalent to saying that 

(3.41) (1 - II) (W A (X, Y) + R(X, Y)N A ) = 0, 

from which (|3~35| follows since W A (X, Y) e TE. □ 
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Note that since 7 -2 P 2 is the projection onto TE one can write (|3.38l) as 
(3.42) ~f 2 W A (X,Y) = -V 2 (R(X,Y)N A ). 

Since both Wa and D x can be expressed in terms of Ba, one obtains the following 
expression for Wa' 

Proposition 3.13. For X, Y e TE one has 

-, 2 W A (X,Y) =(V x Ba)(Y) - (V Y B A )(X) 

- xl 2 )B A {Y) - {Vy1 2 )B a (X) 



+ [g(B A (N B ),X)B B (Y) - g(B A (N B ),Y)B B (X) 



B = l 



As the aim is to express the Codazzi-Mainardi equations in terms of Nambu brack- 
ets, we will introduce maps Ca that is defined in terms of Wa and can be written 
as expressions involving Nambu brackets. 

Definition 3.14. The maps C A ■ C x (E) x • • • x C°°(£) TE are defined as 

(3.43) C A (h, ■ • ■ Jn-2) = Y p eabai "' an ~ 2W A(e a ,e h )(d ai h) ■ ■ ■ (5 a „_ 2 /„_ 2 ) 
for A = 1, . . . ,p and n ^ 3. When n = 2, Ca is defined as 

C A = ^e ab W A (e a ,e b ). 

Proposition 3.15. Let {gi, g 2 }f = {51, 92, fi, ■ . . , fn-2}- Then 

C A (fi,...J n -2y = { 1 - 2 (B A )l,x k } f + ± {x j ,x l } f [r jk (B A )f - (B A )%r^] 

-iE [{< xl } f ( B B)i+nA xl ^ m }f n A^B)l](n B ) k . 

' B = l 

Remark 3.16. In case E is a hypersurface, the expression for C = C\ simplifies to 

C(fu U-2T = {j- 2 B l k , x k } f + ^ {x^x l } f [r jk B? - B*r*,], 
since D X N = 0. 

It follows from Proposition 13.121 that we can reformulate the Codazzi-Mainardi 
equations in terms of Ca '■ 

Theorem 3.17. For all fi,..., f n - 2 e C°°(E) it holds that 

(3.44) 7 2 C A (/ 1 ,...,/„_ 2 ) = (V 2 ))[{x\t{ r } } - {x\x l } f Y^Y{ m y A d i , 
for A = l,...,p, where {51,32}/ = {gi, 92, fi, ■ ■ • , /n- 2 }- 

Proof. As noted previously, one can write the Codazzi-Mainardi equations as 

j 2 W A (X,Y) = -P 2 (R(X,Y)N A ). 
That the above equation holds for all X,Y e TE is equivalent to saying that 

7 2 -^ a6oi "' a "- 2 W A (e a ,e 6 ) = -j-e aba i- a -*V 2 (R(e a ,e b )N A ) 
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for all values of ai, . . . , a n -2 G {1, . . . , n}; furthermore, this is equivalent to 

l 2 C A (h, ■ ■ .,/„_a) = -Y p e aba '- a "- 2 V 2 (R(e a ,e b )N A )(d ai h) ■ ■ ■ (d an _J n . 2 ) 
for all f\, ... , fn-2 £ C' X (S). It is now straightforward to show that 



1 



-s 



abai •••a n — i 



(R(e a , e b )N A )\d ai h) ■ ■ ■ (5 a „_ 2 /„_ 2 ) 



2p 

= ({x k ,T l k: j}f — {x k , x 1 } j f [™f fc TO ^n^, 

which proves the statement. □ 

If M is a space of constant curvature (in which case g(R(X, Y)Z, N A ) = 0), then 
Theorem 13. 171 states that 

(3.45) Ca(/i,...,/„- 2 ) = 

for all /i, ... , /„_ 2 e C X (E). Furthermore, if M = R m , then (|3~44"1) becomes 

(3.46) 7 2 { 1 - 2 (B A y k ,x k } f - j] [{n k A ,x l } f (B B )}](n B ) k = 0. 

B = l 

3.3. Covariant derivatives. Equation (13.25)) tells us that knowing Vx7, for 
X, Y e TE, one can compute VxY through the formula 

VxY = V X Y - ^ E 5(Sa(X),F)7V a , 

^ A = l 

which requires explicit knowledge about the normal vectors. Are there other quanti- 
ties involving V that can be computed solely in terms of the embedding coordinates? 
We will now show that the two derivations 

(3.47) = __j_^ U)5f ' } 

(3.48) V i (u)=g IJ D I {x i )D J (u), 

can be considered as analogues of covariant derivatives on E. Their indices are 
lowered by the ambient metric . Let us start by showing that several standard 
formulas involving covariant derivatives with contracted indices also hold for our 
newly defined derivations. 

Proposition 3.18. For u,v e C X (E) it holds that 

(3.49) Vu = V l {u)d l = D^^D'ix^dt 

(3.50) g(Vu, Vv) = V t (u)V l (v) = Dr^D^v) 

(3.51) A(u) = D t D» = D I D I (u) 

(3.52) |V 2 u| 2 = P.^ftijDjD'fM) = DiD J (u)D J D I (u) 

Proof. The most convenient way of proving the above identities is to work in a 
coordinate system where it 1 ,..., it™ are normal coordinates. In particular, this 
implies that T^ c = 0, which is equivalent to gij(d a x' l )d b 2 c x : > = 0. Let us now prove 
formula (I3.52p for the operators D 1 . 

Let us first note that in normal coordinate one obtains 

|V 2 u| 2 = (V a V b u)(V c y d u)g ac g bd = g ac g bd (d 2 ab u) (d 2 cd u). 
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We now compute 

D 1 D J {u)D J D l {u) = 1 { 7 ~>, x J }, x K }g KI {T l {u, S 1 }, x L }g LJ 
7 (n ij. 

= 2, l —r^ e aS d a {E^{d p u){d^ J )) (dsx^g^d^id^id^ 1 )) (dsx L )g LJ 
g [n ij. 

The terms involving d a dpx J and dcdqX 1 vanish since they appear in combinations 
such as (d a dpX J )(dsx L )g L j which is zero due to the presence of a normal coordinate 
system. Thus, 

D I D J {u)D J D I {u) = 1 e a3 E^g Sg ^e cS gps {dl p u) (8 2 cq u) 
g (n ij. 

= g aq g pc (dl p u)(d 2 cq u) = \v"u\ 2 . 
The other formulas can be proved analogously. □ 

By definition, the curvature tenor of £ arises when one commutes two covariant 
derivatives. In light of Theorem 13.71 one may ask if there is a similar Nambu 
bracket relation which gives rise to the Ricci curvature. A particular example that 
introduces curvature is the following 

(3.53) (V a u)V a V h V b u = (V Q u)V 6 V a V h u - g{1Z(Vu),Vu). 

Since (V a w)V a VbV b u = g(Vu, VAw), it follows from Proposition ^. 18l that one can 
write it as 

(3.54) (V a u)V a V b V b u = V^u^Vp^u) = D I (u)D 1 DjD J '(u), 

and the term in p.53j) involving the Ricci curvature is written in terms of Nambu 
brackets through Theorem 13.71 Using the relation 

(3.55) A(|Vm| 2 ) = 2(V a u)V f, V Q V fc u + 2|V 2 u| 2 , 
and (|3.52[) one obtains 

(V a U )V h V Q V 6U = ]^V{D\V,(u)V\u))-V{D\u)V,V\u) 

where [2? 4 ,2? J ] denotes the commutator with respect to composition of operators. 
Thus, we arrive at the following result: 

Proposition 3.19. Let 1Z be the Ricci curvature o/S and let u e C^E). Then it 
holds that 

V i (u)V i V j T> j (u) = V l {u)V j V 3 V l {u) + [X> i ,D i ](u)Pi'D J '(u) - g(K(Vu), Vu) 
D I (u)D 1 DjD J (u) = D I (u)D J DjD I (u) + {D I ,D J }{u)D I D J '(u) - g(K(Vu),Vu). 

Note that it follows from Theorem 13.71 that the term g(lZ(Vu), Vu) can be written 
in terms of Nambu brackets. If the formulas in Proposition 13.191 are integrated, 
one arrives at expressions whose index structure closely resembles that of equation 
(|3.53p . Namely, by partial integration one obtains 

J (p I {u)D J DjD I {u) + lD I ,D J J(u)D I D J (u)} s fg = J D I (u)D J D I D J {u)^fg, 
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which implies 

(3.56) jD I (u)D I D J Dj(u) y /g = J (d^DjD 1 D J (u) - g(K(Vu), Vu)) y/g. 

Note that since the operators D 1 contain a factor of 7 _1 , the integration is actually 
performed with respect to p, as 7 _1 - v /<7 = p. 

The derivations D 1 and V 1 have indices of the ambient space M; do they exhibit 
any tensorial properties? The object V % {u) transforms as a tensor in the ambient 
space M, i.e. 



v >) = 7 2 (n 1 _i )! { u ^}g/j(y){y < ^} 



i 5y 



,f / }g 7J ( a ;){ a ; fe ,x J } = |^^( U ), 



7 2 (n- 1)! fe fc 1 ' yi ' J dx fc ' 

but this does not hold for the next order derivative V % V 3 (u) due to the second 
derivatives on the embedding functions. One can however "covariantize" this object 
by adding extra terms. 

Proposition 3.20. Define V lJ acting on u e C' X (E) as 

(3.57) V ij (u) = ^(v i V j {u)+V j V i (u)-V u (V i (x j )) S j, 

where T> u (f) = 7 ^^„ 1 _ 1 j ! {/, x I }g IJ {u, x J }. Then V iJ (u) transforms as a tensor in 
M, i.e. 



dy l dy 3 

and for all 1,7 e TE it holds that 



\7 tl (u)X l Y 3 = (V a V b u)X a Y b . 

In particular, this implies that g^V 10 (u) = A (it) and g^g^W^ (u)W 31 (u) = |V 2 u| 2 . 

3.4. Embedded surfaces. Let us now turn to the special case when E is a surface. 
For surfaces, the tensors V, Sa and Ta are themselves maps from TM to TM, 
and Sa coincides with Ta- Moreover, since the second fundamental forms can be 
considered as 2 x 2 matrices, one has the identity 

2 dot W A = ( tr W A ) 2 - tr W\, 
which implies that the scalar curvature can be written as 



R = \(V 2 ) ik (V 2 ) i % jk i+2Y l detW A . 

^ A=l 

Thus, defining the Gaussian curvature K to be one half of the above expression 
(which also coincides with the sectional curvature), one obtains 



p 



(3-58) K='(vy k (vy l R ijkl -^^^ 



' ' A=l 

which in the case when M = M. m becomes 

i p m 

(3.59) K = - J] {jyj&M, 

~ A=li,j=l 
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and by using the normal vectors Z a the expression for K can be written as 



1 _ y Yie jk uE imn i{x l , {x k ,x l }}{x j ,{x m ,x n }} 



r „ Rn , K 8 7 4 (p-l) 
(3-60) 

= ^i^{W,x k },x k }{{x\x l },x l }--{{x^x k },x l }{{x^x k },x 1 } 

To every Riemannian metric on £ one can associate an almost complex structure 
J through the formula 

J(X) = -L £ ac g cb X b e a , 

and since on a two dimensional manifold any almost complex structure is integrable, 
J is a complex structure on S. For X 6 TM one has 

(3.61) V{X) = l —g{X, e a )s ab e b , 

and it follows that one can express the complex structure in terms of V '. 

Theorem 3.21. Defining J M (X) = jV{X) for all X e TM it holds that J M (Y) = 
Sf{Y) for all Y 6 TT,. That is, jV defines a complex structure on TS. 

Let us now turn to the Codazzi-Mainardi equations for surfaces. In this case, the 
map Ca becomes a tangent vector and one can easily see in Proposition 13 . 1 51 that 
the sum in the expression for Ca can be written in a slightly more compact form, 
namely 



C A ={ 1 - 2 (B A )lx k }d l + i {zP,x 1 } [f%(B A )t - {B A )\t k ^ 
+ — J] B b Sa(N b ). 



' B = l 

Thus, for surfaces embedded in R m the Codazzi-Mainardi equations become 



m 1 P 

^ { 7 ~ 2 {x i ,x^}{aP,n k A },x k }d i + : ^ J B b Sa(N b )=0, 

j,k=l 

and in R 3 one has 



7 

j,k=l ' B=l 



(3.62) £ {j- 2 W,x J }W,n k },x k }=0. 
j,k=i 

Let us note that one can rewrite these equations using the following result: 
Proposition 3.22. For M = R m and i = 1, . . . , m it holds that 

m m 

(3.63) 2 {f{x\xi}{x\n k },x k }= £ ^}{^ \ n k } 

j,k — l j,k — l 

for any normal vector N = n l di and any f e C cx ~(Yj). 
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Proof. We start by recalling that for any g e C°°(E) it holds that Ya=x{9i x l }n l = 0, 
since it involves the scalar product g(e a ,N). Moreover, one also has 

£> fe ,n fc } = 2 -e ab (d a x k )(d b n k ) = ^ -e ab (d b (n k d a x k )-n k d 2 ab x k ) 

fc=l k=l P fc=l ^ 

= - 2 wsy = o, 

which implies that Xfcli{ a;fc 'ff" fc } = f° r au 5 6 C X (S). By using the above 
identities together with the Jacobi identity, one obtains 

{f{x\ x*}{x?, n k }, x k } = f{x\ x>}{{x\ n k }, x k } + {x\ n k }{f{x\ X*}, x k } 
= -f{x\x j }{{x k ,x j },n k } -n k {x j ,{f{x i ,x j },x k }} 
= ~f {x l , ^ } { {x k , x> } , n k } + n k { f K , x* } , {x k , x> } } 
= -f{x\x>}{{x k ,xi} 1 n k }-{x k ,xi}{f{x\xi} 1 n k } 
= {f{x\x>}{x\x k },n k }. □ 

Hence, one can rewrite the Codazzi-Mainardi equations for a surface in R 3 as 
(3.64) 2 {^(V 2 y k ,n k }=0, 

j,k=l 

and it is straight-forward to show that 

£ (d c x l ){^- 2 (P 2 y k ,n k } = -e ab V a h bc , 

i.j,k — l 

thus reproducing the classical form of the Codazzi-Mainardi equations. 

Is it possible to verify (|3.64l) directly using only Poisson algebraic manipula- 
tions? It turns out that that the Codazzi-Mainardi equations in K 3 is an identity 
for arbitrary Poisson algebras, if one assumes that a normal vector is given by 

^£ijk{x^,X }di- 

Proposition 3.23. Let {-, •} be an arbitrary Poisson structure on C X (S). Given 
x x ,x 2 ,x z e C rJi '(Y,) it holds that 

2 i £feto {7- 2 K,x^}{^",x fc }, 7 - 1 {^,x"}} = 

j ,k ,1 ,n = l 

for i = 1, 2, 3, where 

2 < 1 2i2 i ( 2 3l2 , ( 3 ll2 
7 = \X ,X } + \X ,X } + {£ ,2 } . 

Proof. Let u,v,w be a cyclic permutation of 1,2, 3. In the following we do not 
sum over repeated indices u, v, w. Denoting by CM 1 the i'th component of the 
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Codazzi-Mainardi equation, one has 

CM" = -{^({x^x"} 2 + {x w ,x u } 2 ), 1 - 1 {x v ,x w }} 

+ {- f - 2 {x u 1 x v }{x v ,x w } ll - 1 {x u ,x v }} + {^ 2 {x u ,x w }{x w ,x v }, 1 ~ 1 {x w ,x u }} 
= -{\- 1 - 2 {x v ,x w } 2 )1 - 1 {x\x w }} +1 - 1 {x^ 1 x v }{ 1 - 1 {x\x w }, 1 - 1 {x u ,x v }} 
+ 1 - 1 {x u ,x w }{ 1 - 1 {x w ,x v } :J - 1 {x w ,x u }} 

= I{ 7 -i{^,^}, 7 - 2 ( 7 2 - {x v ,x w } 2 )} =0. □ 

Let us end by noting that these results generalize to arbitrary hypersurfaces in 
R n+1 . Namely, 

{ 7 - 2 {x\x J }{x J ,n k },x k } f = { 1 - 2 {x\x J }{x J ,x k },n k } f , 

(8 c x l ){ 1 - 2 (P 2 Y k ,n k } = --e ab ^- a ^(V a h bc )(d a Ji)---(d an ^f n ^ 2 ), 
J p 

and 

e klL { 1 - 2 {x\x J }{x J ,x k }, 1 - 1 {x l ,x L }} f = 
for arbitrary x 1 , . . . , x n+1 e C' X (E). 

4. Matrix regularizations 

In physics, "fuzzy spaces" have been used for a long time to regularize quantum 
theories and to model non-commutativity, originating in the study of a quantum 
theory of surfaces (membranes) sweeping out 3-manifolds of vanishing mean curva- 
ture). The main idea was to replace smooth functions on a surface by sequences of 
matrices, approximating the Poisson algebra of functions with increasing accuracy 
as the matrix dimension grows. Since the expressions for geometric quantities de- 
rived in Section [3] uses only the Poisson algebraic structure of the function algebra, 
it is natural to study their matrix analogues in this context. 

Let us start by introducing some notation. Let Nx, iVjj, ... be a strictly increasing 
sequence of positive integers and let T a , for a = 1, 2, . . ., be linear maps from C°°(E) 
to hermitian N a x N a matrices. Moreover, let h : K — > R be a strictly positive 
decreasing function such that limjv^ x Nh(N) converges, and set H a = h(N a ). 
Introduce the operators 

d*(h) = {f,h} 

as well as the matrix operators 

to = i[x,T a (/)], 

and write 

gfi-fk(h) = d h d h ■■■d ik (h) 
t" fk (X) = di 1 d f :-.-d f :(X). 
Let us now define what is meant by a matrix regularization of compact surface. 
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Definition 4.1. Let Ni,N2, ... be a strictly increasing sequence of positive inte- 
gers, let {T a } for a = 1, 2, . . . be linear maps from C X (E,R) to hermitian N a x N a 
matrices and let h(N) be a real-valued strictly positive decreasing function such 
that limjv->oc NH(N) < oo. Furthermore, let w be a symplectic form on E and let 
{-, •} denote the Poisson bracket induced by u. 

If for all integers 1 < I ^ k, {T a } has the following properties for all /, /i, . . . , h e 



C°°(E) 








(4.1) 


Km | 

a— >0O 






(4.2) 


lim | 

Q— >"X 


\T a (fh)-T a (f)T a (h)\\=0, 




(4.3) 


lim 

a— >cc 


d f r h {T a {f))-T a {d^{f)) 


= 


(4.4) 


lim 1nh a Tr (/) = f /w, 





where 1 1 ■ 1 1 denotes the operator norm and h a = h(N a ), then we call the pair (T a , h) 
a C ' -convergent matrix regularization of (E, w). If (T a , h a ) is C k -convergent for all 
k ^ then (T a , h a ) is called a smooth matrix regularization of (E,w). 

In the following, when we speak of a matrix regularization without any reference 
to the degree of convergence, we shall always mean a C 1 -convergent matrix regu- 
larization. 

Remark 4.2. In some cases, a C 1 -convergent matrix regularization is automatically 
a smooth matrix regularization. For instance, if it holds that for any f,he C co (E) 
there exists A/.(f, h) e C rx (E) such that 

^-[T Q (f),T a (h)] = y Ck , a (f,h)T a (A k (f,h)), 

for some Ck, a (f,h) 6 R, then C fc -convergence implies C fc+1 -convergence. The ma- 
trix regularizations for the sphere and the torus in Section 14.21 both fall into this 
category. Hence, they are examples of smooth matrix regularizations. Note that one 
can easily destroy the smoothness of a matrix regularization by slightly deforming 
it, see Example 14.161 

Definition 4.3. A sequence {/ Q } of N a x N a matrices converges to f (or C°- 
converges to f) if 



(4.5) 



lim 

Ct — >GO 



f a -T a (f) 



0. 



Moreover, for any integer k 1, a sequence {f a } of iV Q x N a matrices C k -converges 
to f if in addition 



lim 

a— >x 



f \L)-T a {d h - h {f)) 



0, 



for all 1 ^ I ^ k and /i, • • • , /; e C X (E). If {/ Q } is C^-convergent for all positive k 
then we say that {f a } is a smooth sequence. 

Remark 4.4. If the matrix regularization is C fe -convergent, it is clear that the 
matrix sequence T a (f) is C fe -convergent. It is however easy to construct, even in a 
smooth matrix regularization, C°-convergent sequences that are not C^-convergent; 
see Example 14.151 
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Definition 4.5. A C fc -convergent matrix regularization (T a ,h) is called unital if 
the sequence {1jv q } C fc -converges to the constant function 1. 

Remark 4.6. Although unital matrix regularizations seem natural, and all our ex- 
amples fall into this category, it is easy to construct examples of non-unital matrix 
regularizations. Namely, let (T a ,h) be a matrix regularization and consider the 
map T a defined by 



f a (f) 



( 



Ta(f) 



Then {T a , h) is a matrix regularization which is not unital, since 



lim 

a— >oo 



1. 



T a (l)-t Na + 1 

Proposition 4.7. Let (T a ,h) be a unital matrix regularization. Then 
(4.6) lim 2nN a h a = u. 

Proof. Let us use formula (|4.4I) with / = 1. 

w = lim 27r?i Q Tr Tq,(1) = lim 2nh a Tr T a (X) +liv a — Ijv, 

a— >oo a— »oo L 

lim f 2nh a N a + 2irh a r Tr{T a (l) -livj") = hm 2?rS a 



since 



lim \2irh a Ti(T a (l) -IjvJI < hm 2ith a N a \\T a (l) - t Na 

a— »oo ' a— >x 

due to the fact that the matrix regularization is unital. 



0, 



□ 



Proposition 4.8. Let (T a ,h a ) be a C k -convergent matrix regularization and as- 
sume that f a and h a C k -converge to f,he C X (S) respectively. Then it holds that 
a fa + bh a C k -converges to af + bh, for any a, b e R, and / a /i Q C k -converges to fh. 
Furthermore, it holds that 



(4.7) 
(4.8) 



lim 

a— »oc 



/a 



lim ||T a (/)|| 

a— >QC 

lim 27r?i a Tr (f a h a ) = fhu). 

K-VOO J s 



Proof. The fact that a/ + 6/i C^-converges to a/ + 6/i follows directly from linearity 
of the maps T a . To prove (|4.7I) one uses the reverse triangle inequality to deduce 



lim 

a— >rc 



H/a||-||T a (/)|| 



^ lim 

a— »oc 



f a -T a (f) 
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since f a is assumed to converge to /. Let us continue by proving that f a h a C°- 
converges to fh, i.e. 



lim 

a— >oo 



faha. ~ T a (fh) 
lim ( 

a— >oo \ 



lim 

a — »O0 



f a h a - f a T a (h) + f a T a (h) - T a (fh) 



sC lim 

a— >oo 



/a 






+ 


f a T a (h) - T a (f)T a (h) + T a {f)T a (h) 


/a 




- r a (ft) 


+ 


L-T a (f) 


\\T a (h)\\ + \\T a (f)T a (h) 



T a (fh)\\) 
Ta(fh)\\) 



0, 



since both {f a } and {h a } are C°-convergent sequences and is bounded by 

(|4.7[) . Using the face that f a h a C°-converges to fg, it is easy to prove (|4.8[) by 
computing 

lim 2nh a Tr f a h a = lim 27rS a Tr (f a h a - T a (fh) + T a (fh)\ 

a— >X' a— »oo 

= lim 2irh a TrT a (fh)) = fhui. 

Finally, we proceed by induction to show that f a h a C fc -converges to fh. Thus, 
assume that, for some ^ I < k, u a v a C'-converges to uv whenever u a and v a 
C'-converges to u and v respectively. Since 

da {f a h a ) = (d^f a )h a + f a d^h a 

we can use the induction hypothesis (together with the assumption that f a ,h a 

C fc>i -converges) to conclude that (f a h a ) C'-converges, which implies that f a h a 
C l+1 -converges. Hence, it follows that f a h a C fe -converges to fh. □ 

The above result allows one to easily construct sequences of matrices converging to 
any sum of products of functions and Poisson brackets. Namely, simply substitute 
for every factor in every term of the sum, a sequence converging to that function, 
where Poisson brackets of functions may be replaced by commutators of matrices. 
Proposition 14.81 then guarantees that the matrix sequence obtained in this way 
converges to the sum of the products of the corresponding functions, as long as the 
appropriate level of convergence is assumed. 

Proposition 4.9. Let (T Q , h) be a matrix regularization and let {/«} be a sequence 
converging to f . Then lim a ^ x ||/ a || =0 if and only if f = 0. 



Proof. From Proposition 14.81 it follows directly that if f a converges to then 

lim ||/ Q || = lim ||T Q (0)|| =0. 

a— >oo a— »cc 



Now, assume that lun^rx, ||/ Q || = 0. Then it holds that 
j r- , _ i;,„ •).--. tv r- i;... •>--/, v II r- 

from which we conclude that / = 



f 2 uj = lim 2-Kha^tfi ^ lim 2irh a N a \\fl\\ ^ lim 2nh a Na\\f a 

Oi— >GC a—>c/j CK — >OC 



0. 



□ 



Proposition 4.10. Let (T a ,h) be a matrix regularization and assume that {f a 
C k -converges to f. Then {f a } C k -converges to f. 
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Proof. Due to the fact that 



1 1 -At 1 1 one sees that 



lim 

ol — >0C 



lim 

a— >x 



lim 

a—*cc 



d f r f "(f a )-T a (d f i- f «(f)) 



since {f a } C fe -converges to /. 



□ 



Proposition 4.11. Let (T a ,h) be a unital matrix regularization and assume that 
f is a nowhere vanishing function and that {f a } C k -converges to f. If f^ 1 exists 
and H/a 1 )! is uniformly bounded for all a, then {Z^ 1 } C k -converges to 1//. 

Proof. Let us first show that f^ 1 C°-converges to 1//; one calculates 



lim 

a— >od 



f-^-T a {l/f) 



^ lim 

a— i>cc 



= lim 

a— >cc 

^ lim 
= 0, 



r 1 
r 1 



fa 1 lN a -f a T a (l/f) 

f a T a (l/f)+T a (l)-T a (l) 



L N a 



■T a (l) 



/«r a (i//)-r a (i)||) 



since the matrix regularization is unital and H/q 1 )! is assumed to be uniformly 
bounded. Let us now proceed by induction and assume that f^ 1 is C'-convergent 
(0 s£ I < k). For arbitrary h e C°°(£) it holds that 

[f-\T a {hJ\ = -f^tfcMh)]?- 1 , 

and since f a is C fc -convergent, the above sequence is C'-convergent by Proposition 



14.81 which implies that f a is C + -convergent. Hence, it follows by induction that 
/qT 1 is C fc -convergent. □ 

4.1. Discrete curvature and the Gauss-Bonnet theorem. Let us now con- 
sider a surface E embedded in M via the embedding coordinates x ,x m , with 
a symplectic form 

u) = p(u 1 ,u 2 )du 1 a du 2 , 

inducing the Poisson bracket {/, h] = j;£ ab (d a f)(dbh), and let (T a ,h a ) be a matrix 
regularization of Furthermore, we let {7 a } be a C 2 -convergent sequence 

converging to 7 = -Jg/p (and we assume that {7a 1 } exists and converges to I/7), 
and we set X l a = T a (x l ) as well as N Aa = T a (n l A ) for i = l,...,m. Moreover, 
given the metric and the Christoffel symbols t l - k of M, we let {GV,-, Q } and 

{ij/c a ) denote sequences converging to and T* fc respectively. To avoid excess of 
notation, we shall often suppress the index a whenever all matrices are considered 
at a fixed (but arbitrary) a. 

Since most formulas in Section [3] are expressed in terms of the tensors Pj and 
(SaYj (in the case of surfaces), we introduce their matrix analogues 



1 



ih 
1 



j'j 



23 



as well as their squares 

(V 2 )) = (n)^ and {& a )) = {S a WSa$, 
and corresponding trace 

m m 

CvV^^CP 2 )} and trSi = £(Si)*. 

i=l i=l 

(The ordinary trace of a matrix X will be denoted by Tr X.) From Proposition ^. Sl it 
follows that one can easily construct matrix sequences converging to the geometric 
objects in Section [31 as long as the appropriate type of convergence is assumed. 
Let us illustrate this by investigating matrix sequences related to the curvature of 
S and the Gauss-Bonnet theorem. 

Definition 4.12. Let (T a ,h) be a matrix regularization of (£,w), let K be the 
Gaussian curvature of £ and let x be the Euler characteristic of A Discrete 
Curvature of £ is a matrix sequence {K\, K2, . . .} converging to K, and a Dis- 
crete Euler Characteristic ofT, is a sequence {xi, Y2, X3j • ■ •} such that lim Xa = X- 

a— »oo 

From the classical Gauss-Bonnet theorem, it is immediate to derive a discrete ana- 
logue for matrix regularizations. 

Theorem 4.13. Let (T a , h) be a matrix regularization o/(£, uj), and let {Ki, K 2l ■ ■ ■} 
be a discrete curvature of S. Then the sequence X11X2, ■ • ■ defined by 

(4.9) Xa = h a Tr \%K a , 
is a discrete Euler characteristic o/S. 

Proof. To prove the statement, we compute lim Q _> x Xa and show that it is equal 
to x(S). Thus 

lim Xa = lim ^-2nh a Tr \%K a , 

a — »co a— >x Z7T L 

and by using Proposition I4.8I we can write 

lim Xa = — f K^-uj = — [" K^-pdudv = — f K-Jgdudv = x(£), 

a->cc 27T J s p 27T J s p 2lT J s 

where the last equality is the classical Gauss-Bonnet theorem. □ 

Theorem 4.14. Let (T a ,H) be a unital matrix regularization of (£,w) and let 
Rijki, for each i,j,k,l = l,...,m, be a sequence converging to the component of 
the curvature tensor of M . Then the sequence K defined by 

Z A = l 

is a discrete curvature o/E. Thus, a discrete Euler characteristic is given by 

(4.10) x = STr (r 3 (r 2 r k (V 2 y l R l3kl ) - % j] Tr \r^S 2 A . 

Z A=l 1 J 

Proof. By using the way of constructing matrix sequences given through Proposi- 
tion [4J2 the result follows immediately from Theorem 13.71 □ 
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In the case M = R m it follows from the results in Section 13.41 that when (T a , K) is 
a C 2 -convergent matrix regularization, then the sequence 



(4.11) 



(7l ) 2 [[xi,X*\,XZ\ [[XI, X l a ], Xi]7 Q 2 

1 (7l) 2 [[XL x k a ],x l a ] [[xi, X*],X l a ]% 



4 

converges to the Gaussian curvature of E. 
4.2. Two simple examples. 

4.2.1. The round fuzzy sphere. For the sphere embedded in R 3 as 
(4.12) x = (x l ,x 2 , x 3 ) = (cos tp sin 9, sin ip sin 9, cos 9) 

with the induced metric 

^ <**>=(£ sin 2 . 

it is well known that one can construct a matrix regularization from representations 
of su(2). Namely, let Si, S2, S3 be hermitian N x N matrices such that [S J , S k ] = 
iei k iS\ (S 1 ) 2 + (S 2 ) 2 + (S 3 ) 2 = (iV 2 - l)/4, and define 

(4.14) X 1 = . 2 S\ 

1 ' V^ 2 -1 

Then there exists a map (which can be defined through expansion in spherical 
harmonics) such that TW(x') = X 1 and (T^ N \h = 2/^jN 2 - 1) is a unital matrix 
regularization of (S 2 ,^/gd9 a d<y9) GH82 . A unit normal of the sphere in R 3 is 
given by N e TR 3 with N = x i di, which gives N l = X 1 , and one can compute the 
discrete curvature as 

(4.15) ^ N = ~Jf 2 Tr[X i ,X i ] 2 =lit 

i<j=l 

which gives the discrete Euler characteristic 

27V 



(4.16) XN = hTrK N = hN- 

\JN Z — 1 

converging to 2 as 2V — > 00. 

4.2.2. The fuzzy Clifford torus. The Clifford torus in S 3 can be regarded as em- 
bedded in R 4 through 

x = (x 1 , x 2 , x 3 , x 4 ) = — = (cos <pi, sin <pi, cos <p% , sin (^2) , 
V2 



with the induced metric 



Jab) - 2 (o ! 



and two orthonormal vectors, normal to the tangent plane of the surface in TR , 
can be written as 

N+ = x 1 d 1 + x 2 d 2 + x 3 d 3 + x 4 d 4 . 
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To construct a matrix regularization for the Clifford torus, one considers the N x N 
matrices g and h with non-zero elements 

gkk=u> k ~ 1 for k = 1,...,N 
hk,k+i = 1 for k = 1, . . . , N - 1 

where to = exp(i29) and 9 = ir/N. These matrices satisfy the relation hg = ujgh. 
The map is then defined on the Fourier modes 



as 

and the pair (T^ N \ h = sin 9) is a unital matrix regularization of the Clifford torus 
with respect to ^Jgdipi a d<f2 |FFZ89| |Hop 89 . Thus, using this map one finds that 

' (5 f +5) 



X 1 = 


TV) 


= ^ T (cosv5i) 


X 2 = 


T{x 2 ) 


= ^ T ( sin< ^ 1 ) 


X 3 = 


T(x 3 ) 


= -^ T (cosv5 2 ) 


x 4 = 


T(x 4 ) 


= ^ T ( sin( ^2) 


which implies that N\ 


= x\ 


Nl = X 2 , Nl 



2V2 

: 2 , Nl = +} 

straightforward computation one obtains 
and therefore 



2V2 



and since [X\X 2 ] = [X 3 ,X 4 ] = it follows that 

2^2 t [X\Ni][X^,Nt] = ^- 2 t t Xi > Xj ? = - 1 - 

i,j = l j,J = l 

This implies that the discrete curvature vanishes, i.e. 

i,j=l »,J=1 

which immediately gives = 0. 

The following two examples will show that even in the smooth matrix regular- 
ization of the torus it is easy to find sequences that are not smooth, and that the 
regularization can be deformed into a non-smooth matrix regularization. 
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Example 4.15. Let (T a ,h a ) be the matrix regularization of the Clifford torus as 
in Section ^. 2. 2\ For each N, define the matrix 

= diag(/i s ,O,...,O), 

for some fixed < s ^ 1 . Clearly, it holds that 



lim 

a— >oo 



T a (0) 



lim 

a— »oc 



i.e. 9 C -converges to 0. Let us show that 9 does not C 1 -converge to 0. If 
C 1 -converges to 0, then it must hold that 



lim 

a— >oo 



_[0,T Q (/)]-T Q ({O,/}) 



lim 

a— »cc 



for all f 6 C X (S). For i? = 2V2T (Ar) (a: 3 ) 



one computes the eigenvalues 



ofA = ±[9,H] to be 

Ai = i^h 8 - 1 A 2 = -i\j2h a ~ l A 3 = • • • = Aat = 0. 

Hence, the norm of A does not tend to 0, which implies that 9 is not C 1 -convergent. 

Example 4.16. Let (T a ,h a ) be the matrix regularization of the Clifford torus as 
in Section \4-.2.S\ For each N, define the matrix 

= diag(/i s ,O,...,O), 

for some fixed 1 < s ^ 2. Let us now deform the fuzzy torus to obtain a C 1 - 
convergent matrix regularization that is not C 2 -convergent. Defining 

S a (f)=T a (f) + n(f)9, 

where fi : C°°(£) — » R is an arbitrary linear functional, one can readily check that 
(S a ,H a ) is a C 1 -convergent matrix regularization of the Clifford torus. Let us now 
prove that (S a ,h a ) is not a C 2 -convergent matrix regularization, and let us for 
definiteness choose fi to be the evaluation map at ifi = if2 = 0. 
In a C 2 -convergent matrix regularization it holds that 



lim 

a— >oo 



h 2 



[[S a (u), S a (v)], S a (w)] - S a ({{u, v}, w}) 



for all u,v,w e C°°(S). Choosing u = 2^/2 cos if2 and v = w = 2s[2 sin if2 gives 
S a (u) = tf + h + 2^/29, S a (v) = i{tf - h) and {u, v} = 0. Thus 



lim 

a— >rc 



^[[S a (u),S a (v)],S a (w)] - S a ({{u,v},w}) 

2V2 



[[§, - h)],i(hl - h)] = lim 2V2(2 + V6)/i s 



which does not converge to 0. Hence, (S a ,h a ) is a C 1 -convergent, but not C 2 - 
convergent, matrix regularization of the Clifford torus. 
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4.3. Axially symmetric surfaces in M 3 . Recall the classical description of gen- 
eral axially symmetric surfaces: 

(4.17) x = (f(u) cosv, f(u) smv, h(u)) 

H = ~~rur^=^77=^ i h '^ cosv > h '( u "> sin v ' -/'("))' 

which implies 

7' 2 +fr' 2 0\ . ±1 fh'f"-h"f 



{9ab) = I q /2 J {hab) 



where /i a h are the components of the second fundamental form. The Euler charac- 
teristic can be computed as 

1 f P+ h'(h'f" -h"f) f> 

which is equal to zero for tori (due to periodicity) and equal to +2 for spherical 
surfaces (/'("+) = +oo if u = h). 

While a general procedure for constructing matrix analogues of surfaces embed- 
ded in R 3 was obtained in [ABH+ 09bj I ABH+ 09a) (cp. also |Arn08bj ). let us restrict 
now to h(u) = u = z, hence describe the axially symmetric surface E as a level set, 
C = 0, of 

(4.19) C {x)= l -{x 2 +y 2 -f 2 {z)) 1 

to carry out the construction in detail, and make the resulting formulas explicit. 
Defining 

(4.20) {F{x), G(x)} K 3 = VC ■ (VF x VG) , 
one has 

(4.21) {x, y} = -ff'{z), {y, z) = as, {z, x} = y, 
respectively 

(4.22) [X, Y] = ihff'(Z), [Y, Z] = ihX, [Z, X] = ifrY 

for the "quantized" ( "non-commutative" ) surface. In terms of the parametrization 
given in (|4.17[) . the above Poisson bracket is equivalent to 

(4.23) {F(u,v),G(u,v)} = e ab {d a F){d b G) 

where d\ = d v and 82 = d u . By finding matrices of increasing dimension satisfying 
(|4.22l) , one can construct a map T a having the properties (|4.2[) and (|4.3[) of a matrix 
regularization restricted to polynomial functions in x, y, z (cp. Ar n08a| ). 

For the round 2-sphere, f(z) = 1 — z 2 , (I4.22[) gives the Lie algebra su{2), and its 
celebrated irreducible representations satisfy 

(4.24) X 2 + Y 2 + Z 2 = 1 if h 



V^ 2 -1' 
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When / is arbitrary, one can still find finite dimensional representations of (|4.22|) 
as follows: rewrite (|4.22p as 

(4.25) [Z, W] = hW 

(4.26) [W, W*] = -2hff'(Z) 

implying that Zi—Zj = h whenever Wij ^ and Z diagonal. Assuming W = X +iY 
with non-zero matrix elements Wk,k+i = w k for k = 1, . . . , N — 1, one thus obtains 
(with wq = wn =0) 

Z kk = |(JV + l-2fc) 

wt - w 2 k _, = -2hff'(h(N + 1 - 2fc)/2) = Q k , 
which implies that 

k 

< = 2 Qi 

1=1 

and the only non-trivial problem is to find the analogue of (|4.24j) . To this end, 
define 

(4.27) f 2 = X 2 + Y 2 = i (WW^ + W^W) , 

with W given as above. As Z has pairwise different eigenvalues, the diagonal matrix 
given in (|4.27p can be thought of as a function of Z; hence as f 2 (Z). It then trivially 
holds that 

(4.28) C = X 2 +Y 2 — f 2 (Z) = 0, 

for the representation defined above. The quantization of h comes through the 
requirement that f 2 should correspond to / 2 . While for the round 2-sphere f 2 
equals f 2 , provided h is chosen as in (|4.24[) . it is easy to see that in general they 
can not coincide, as 

[X 2 + Y 2 -f(Z) 2 , W] = [(WW^ + W^W)/2 - f(Z) 2 , W] 

= \W{W\ W] + \\W\ W]W - f(Z)[f(Z), W] - [f(Z), W]f(Z) 

= ■■■ = f(z)(hf(z)w - [f(Z),w]) + (hf(z)w - [f(z),w])f(Z) 

with off-diagonal elements 

(f(z k ) + /(z fc _i)) (hf'(z k ) - (f(z k ) - /(%_i))) 

that are in general non-zero (hence X 2 + Y 2 + f 2 (Z) is usually not even a Casimir, 
except in leading order). 
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How it does work is perhaps best illustrated by a non-trivial example, f(z) 
l-z 4 : 

h 4 / 

(4.29) w\ = —{(N + ifk - 3(N + l) 2 k(k + 1) + 

2(N + l)k(k + l)(2k + 1) - 2k 2 (k + l) 2 ) 

ft = \H + = ^ (V + l) 3 (2fc - 1) - 6(N + l) 2 k 2 

+ A(N + l)k(2k 2 + 1) - 4fc 2 (fc 2 + 1)) 

(note that w 2 = w 2 N = is explicit in (|4.29l0 so that 

(iV + 1) 4 (N + lf 



(4.30) (X 2 + Y 2 + Z 4 ) kk = h 4 



16 



+ k(N + l)-fc 2 



Expressing the last two terms via Z 2 (note that the cancellation of k 3 and k 4 terms 
shows the absence of Z z and higher corrections) one finds 



ft ^ ((TV + l) 2 - 4(iV + 1) + 4)l 



16 



which equals 1 if ft is chosen as 2/\jN 2 — 1. Note that this is the same expression 
for ft then for the round sphere, f 2 = 1 — z 2 (cp. (|4.24p ). 

A more elegant way to derive the quantum Casimir (cp. also |Roc91[ [GPS09 ) 

(4.31) Q = X 2 + Y 2 + Z 4 + h 2 Z 2 

is to calculate 

[X 2 + Y 2 + Z 4 , W] = [(WW* + W*W)/2 + Z 4 , W] 

= ... = H 2 [W,Z 2 ], 

which determines the terms proportional to ft in the Casimir. 
Due to the general formula 

(4-32) K = --L £]kl e lpq (^ )- 2 [X\ [X k , X 1 ]] [X\ [X^X^T 2 

one obtains, for the axially symmetric surfaces discussed above, 

(4.33) K = T 2 {iff'f(Z) + ±[W, ff'{Z)]W* + ^[W, ff'(Z)])r 2 
with 

(4.34) 7 2 = i (WW* + W^W) + {fff{Z) = f{Z) 2 (f'(Z) 2 + 0(H), 
giving 

(4.35) k = -{f(zf + t)- 2 f(z)-\nz) + 0(h) 

and for f(z) 2 = 1 — z 4 one has 

(4.36) K = (4Z 6 + 1-Z 4 ) ~ 2 (6Z 2 - 2Z 6 ) + 0(h) 

(4.37) 7 2 =t-Z 4 +AZ 6 + 0(h). 
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Note that (cp. P~55jt ) 



H for arbitrary /, and that (due to the axial 



symmetry) K and 7 are diagonal matrices, so that 



in this case simply being a Riemann sum approximation of ^ K^Jg, indeed converges 
to 2, the Euler characteristic of spherical surfaces. 



4.4. A bound on the eigenvalues of the matrix Laplacian. As we have shown, 
many of the objects in differential geometry can be expressed in terms of Nambu 
brackets. Let us now illustrate, in the case of surfaces, that some of the techniques 
used to prove classical theorems can be implemented for matrix regularizations. 
In particular, let us prove that a lower bound on the discrete Gaussian curvature 
induces a lower bound for the eigenvalues of the discrete Laplacian. For simplicity, 
we shall consider the case when M = M m and, in the following, all repeated indices 
are assumed to be summed over the range 1, . . . , m. 

Let us start by introducing the matrix analogue of the operator D l : 

D a (X) = -±-j- 1 [X,X a ]. 

These operators obey a rule of "partial integration" , namely 
(4.38) Tr (%D a (X)Y) = - Tr (%D a (Y)X) , 

which is in analogy with the fact that 

f (rfD i (f)h)x> = - f {lD\h)f)u. 
In view of Proposition 13. 181 it is natural to make the following definition: 

Definition 4.17. Let (T a ,h a ) be a matrix regularization of (S, u). The Discrete 
Laplacian on £ is a sequence {A Q } of linear maps defined as 

A a (X) = DiDi(X) = -L^[^[X,X a ],X a ], 

where A is a N a x N a matrix. An eigenmatrix sequence of A a is a convergent 
sequence {u a } such that A a (u a ) = X a u a for all a and lim A„ = A. 

a— >oo 

Proposition 4.18. A C 2 -convergent eigenmatrix sequence of A a converges to an 
eigenfunction of A with eigenvalue A = lim A a . 

a— »co 

Proof. Given the assumption that u a is a C 2 -convergent matrix sequence converg- 
ing to u, we want to prove that Au — Xu = 0. By Proposition ^. 101 this is equivalent 
to proving that lim Q ,_ >a0 ||T Q (Au — Au)|| = 0. One obtains 



lim \\T a (Au-Xu) 

a— >cc 



lim 

a— >GO 



T a (Au) - A a U a + A a U a - XT a (u) + Xii a - Xil a 



lim 

a— >x 

= lim 

a— »cc 



Ta{Au)- Aaila + \X\\\-T a {u) + U a \\ + 



^ lim 

Ol — >GO 



+ \X- X a \ \\u a \ 



0, 



since A a u a — X a ii a = and A Q converges to A. 



□ 
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The way curvature is introduced in the classical proof of the bound on the eigen- 
values, is through the commutation of covariant derivatives. Let us state the cor- 
responding result for matrix regularizations. 

Proposition 4.19. Let (T Q , h a ) be a C 2 -convergent matrix regularization o/(£, uj). 
//{«„) is a C 3 -convergent matrix sequence then 



lim 

a— >x 



D^u a )D l a bi a bi a {u a ) - D l a (u a )DiDiD l a (u a ) 

- lKDU(u a )KDi(u a ) + K a D l a (u a )& a (u a ) 



0, 



where [•, •] denotes the commutator with respect to composition of maps. 

Proof. The result follows immediately from Proposition 13.191 and Proposition 14.81 
Note that in the case of surfaces it holds that 1Z a b = Kg a b, where K is the Gaussian 
curvature of S. □ 

A useful corollary is the following: 

Proposition 4.20. Let (T a , H a ) be a C 2 -convergent matrix regularization o/(S, w). 
If {ua} is a C 2 -convergent matrix sequence then 

lim S a Tr (%bl l bibi(u a )bi x (u a )) = 

a— >oo \ / 

lim h a Tr (%bibibi{u a )Di(u a ) - %K a fy(ii a )t>i(ii a )) 

a— >cc \ / 

Proof. It follows from Proposition 14.191 that for a C 3 -convergent sequence u a it 
holds that 

lim ft a Tr (%DiDiDi(u a )Di(u a ) -7a^^j,(6a)^j,(«a) 

a— >oo V 

-7aI^,^](^)^^i(«a)+7a-K'a^(«a)^(«a)) = 0. 

Due to the appearance of a trace, the above holds even for C 2 -convergent sequences, 
since e.g. 

n a Tr%bibibi(u a )bi(u a ) = -h a T±%bibi(u a )bibi(u a ), 

and the latter expression only requires C 2 -convergence. Thus, one obtains 

lim^Tr (%bibibi(u a )bi(u a )) = lim h a Tr (^ a blblb l a (u a )b l a (u a ) 

+ %{b^bi\{ Ua )b\bi(u a ) - %K a K(u a )bi(u a )) 

= iim n a Tv (%bib t a bi(u a )b i a (ua) - %k a b i a ( Ua )b i a (u a )\, 

by using equation (|4.38|) . □ 

Proposition 4.21. Let (T a ,h a ) be a matrix regularization o/(E,cj). If {u a } is a 
C 2 - convergent matrix sequence then 

lim h a Tr (bibi(u a )bibi(u a )) ^ I lim h a Tr (A a (u a )) 2 . 

a— >cc \ J Z a— >oo v ' 
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Proof. By using the fact that |V 2 m| 2 ^ i(Au) 2 (for 2-dimensional manifolds) one 
obtains 

Urn h a Tr (# a ^(u Q )^Di(u Q )) = ^ J |V 2 u| 2 ^ ^ ^ J (Ati) 2 u; 

1-2 

= lim -S a Tr (A Q (u Q )) , 

a— »oe Z 

since M a is assumed to C 2 -converge to u. □ 

Theorem 4.22. Let (T a , h a ) be a C 2 -convergent matrix regularization of (E,u;) and 
fei j« a } be a C 2 -convergent eigenmatrix sequence of A Q iirat/i eigenvalues {— A Q }. 
If K a J; k1at q /or some keM and aZZ a > ao, £/ien lim A Q >- 2k. 

a— >oc 

Proof. Let {u a } be a hermitian eigenmatrix sequence of A a with eigenvalues {— A Q }. 
First, one rewrites 

Tr 7Q A Q (u Q ) 2 = Tr (%D l a D\(u a )DlDl{u a )) 
(4.39) „. „ 

= — A Q Tr (fi Q %£>J,I£(fi a )) = A Q Tr (*f a D l a {u a )D' a {u a )) . 

Then, one makes use of Proposition 14.201 to write 

lim/i Q Tr7 Q A Q (u Q ) 2 = - lim h a Tr {%D l a b' a D^{u a )fD l a {u a )) 

a— ¥00 a— »ao v y 

= lim na^f-7a^a^(«a)^U«a)+7a^a^a(«a)^(«a)) 
a— »oo V / 

= lim ^Trf 7a ^^(u a )^^(uo)+7o^a^(w a )^(u a )y 

a— >oo V / 

Using the assumption that K a ^ ni together with Proposition 14 . 2 1 1 one obtains 

"\2 !• fc rp I 1 - A /„», \2 



lim fi, a Tr7 Q A Q (u ct ) > hm S Q Tr [ -%A a (u a ) + K%D l a (u a )D l a (u a ) 

ot — >oc a— »cc \ z 



Jim Qa q + K^h a Tr ( 7Q £>^(u a ).D^(u Q )) , 



where (|4.39p has been used. One can now compare the above inequality with (|4.39[) 
to obtain 

i(A - 2k) lim /i Q Tr (%D l a (u a )D l a (u a )) ^ 0. 

2 a— >oo v * 



Since 



lim h a Tr (%D l a (u a )D l a {u a )) = ±- f 7 |V U | 2 w ^ 0, 

a— >0O Z7T J s 



due to the fact that 7 is a positive function, it follows that A ^ 2k. □ 

Although the above proof depends on the fact that the matrix regularization is 
associated to a surface (and therefore, the results of differential geometry can be 
employed), we believe that, under suitable conditions on the matrix algebra, there 
exists a proof that is independent of this correspondence. 
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